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ABSTRACT 

Context. Using the multi-integral-field spectrograph GIRAFFE at VLT, we previsouly derived the stellar-mass Tully-Fisher Relation 
(smTFR) at z ~ 0.6 for a representative sample of 63 emission-line galaxies. We found that the distant relation is systematically offset 
by roughly a factor of two toward lower masses from the local relation. 

Aims. We extend the study of the evolution of the TFR by establishing the first distant baryonic TFR in a CDFS subsample of 35 
galaxies. We also investigate the underlying cause of the large scatter observed in these distant relations. 

Methods. To derive gas masses in distant galaxies, we estimate a gas radius and invert the Schmidt-Kennicutt law between star 
formation rate and gas surface densities. We consider the influence of velocity dispersion on the scatter of the relation, using the 
kinematic tracer 5 suggested by Kassin and collaborators. 

Results. We find that gas extends farther out than the UV light from young stars, a median of ~30%. We present the first baryonic 
TFR (bTFR) ever established at intermediate redshift and show that, within an uncertainty of ±0.08 dex, the zeropoint of the bTFR 
does not appear to evolve between z~0.6 and z=0. On the other hand, we confirm that the difference between the local and distant 
smTFR is significant, even considering random and systematic uncertainties, and that accounting for velocity dispersion leads to a 
significant decrease in the scatter of the distant relation. 

Conclusions. The absence of evolution in the bTFR over the past 6 Gyr implies that no external gas accretion is required for distant 
rotating disks to sustain star formation until z=0 and convert most of their gas into stars. Finally, we confirm that the larger scatter 
found in the distant smTFR, and hence in the bTFR, is caused entirely by major mergers. This scatter results from a transfer of energy 
from bulk motions in the progenitors, to random motions in the remnants, generated by shocks during the merging. Shocks occurring 
during these events naturally explain the large extent of ionized gas found out to the UV radius in z~0.6 galaxies. All the results 
presented in this paper support the "spiral rebuilding scenario" of Hammer and collaborators, i.e., that a large fraction of local spiral 
disks have been reprocessed during major mergers in the past 8 Gyr. 

Key words. Galaxies: evolution; Galaxies: kinematics and dynamics; Galaxies: high-redshifts; galaxies: general; galaxies: interac- 
tions; galaxies: spiral. 



1. Introduction 

The stellar-mass Tully-Fisher Relation (smTFR) has received 
increased attention over the past decade , both at low (e.g. , 
iBell & de Jongl |2Q0lb iPizagno et al I 120071: [Mever et alJ [2008) 
and high redshift dConselice et al.l 120051: iFlores et all 2006; 
Atkinson et ail l2007t iKassin et alj [20071 IPuech et al.l [2008 



Cresci et al. I2009I) . However, the smTFR does not allow us to 



homogene o usly c haracterize all galaxy morphologies together. 
iMcG augh (2000) showed that dwarf galaxies, which have a 
larger gas content that more massive galaxies, fall downward the 
smTFR. Considering the baryonic (i.e., stellar plus gas) content 
of galaxies, he showed that a ll galaxies, includin g dwarves, de- 
fine a common TFR (see also Beg um et alj |2008). which by then 
was dubbed the baryonic TFR (bTFR). Because the bT FR ap- 
pears to hold over five decades in mass dMcGaughll2000l) . it can 
therefore be considered as being somehow more "fundamental" 
than the smTFR. 

So far, the bTF R has b e en st u died only in the lo- 
cal Universe ( see l^rheijenl l2Q01j [gell & de Jongl 120011: 
Goruvich et all l2004t iMcGaughl 120041 l2005t lGehaetal.1 



20061: iNoordermeer&Verheiienl 120071: iDe Riicke et al.i 
Avila-Reese et al.l 120081: iMever et al.l l2008t IStark et all 



2007 



2009; 



Trachternach et al.l l2009). This was primarily because HI obser- 



vations of galaxies are practically limited to z<0.3 dLah et all 
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120071) . To estimate the gas content of more distant galaxies, 
indirect methods were develo ped, such as the inversion of the 
Schmidt-Kennicutt (SK) law dKennicutll [T989). which relates 
the gas surface and star formation r ate (SFR) densities (e.g., 
lErb et al.|p2006l: iMannucci et al.ll2009l) . Another difficulty in es- 
tablishing a reliable bTFR (or even a smTFR) at high z is that it is 
not always very clear whether all distant samples are truly repre- 
sentative of the luminosity or mass function of galaxies at those 
epo chs. At intermedia te redshifts (i.e., z<l). lFlores et al] (T2006) 
and I Yang et al.l d2008l) assembled a representative sample of 63 
z~0.6 emission-line galaxies observed by 3D spectroscopy for a 
project called IMAGES. The kinematics of the entire galaxy sur- 
face is probed by 3D spectroscopy, which allows us to classify 
them as a function of their relaxation state. Rotating Disks (RDs) 
are galaxies that exhibit regular rotation well-aligned along the 
morphological axis. Perturbed rotators (PRs) have large-scale 
rotation with a local perturbation in the velocity dispersion map 
that cannot be accounted for by rotation, while galaxies with 
complex kinematics (CKs) do not exhibit large-scale rotation, or 
have a stron g misalignment be tween the dynamical and morpho- 
logical axes (Yang et al. 2008). Using this new kinematical clas- 
sification, Flor es et all ((2006) showed that the larger dispersion 
of the distant TFR was due entirely to galaxies with non-relaxed 
kin ematics, pro b ably a ssociated with mergers, as lat er confirmed 
by IPuech et al l (120081) (hereafter P08) as well as IKassin et all 
(20071) and lCovington et all (12009). Restricting the smTFR to ro- 
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tating disks (RDs), for which the TFR can be confidently estab- 
lished, for the first time P08 detected an evolution in zeropoint of 
the smTFR by dex between z~0.6 and z=0. However, 

z~0.6 galaxies have a larger gas fraction than their local coun- 
terparts (i.e., <f g as >~30% at z~0.6-0. 8), as derived from th e 
evolution of the gas-metallicity relation (Rodrigues et al. 2008). 
This led P08 to suggest that the bTFR, as opposed to the smTFR, 
might not be evolving. 

In this paper, we take advantage of the representative 
IMAGES sample of intermediate-mass galaxies at z~0.6 to es- 
tablish the first bTFR at high redshift, using the SK inversion 
method. The paper is organized as follows: In Sect. 2, we present 
the sample and the data used in this paper; In Sect. 3, we extend 
the previous analysis of the smTFR by P08, obtaining new in- 
sight about its shift in zeropoint and its scatter; In Sect. 4, we 
estimate the gas content of intermediate-redshift galaxies and 
compare it to the stellar light from young stars; In Sect. 5, we 
establish the first bTFR at z~0.6, which is discussed in Sect. 
6. Throughout, we adopt Hq = 70 km/s/Mpc, Qm = 0.3, and 
Qa = 0.7, and the AB magnitude system. 

2. Data & sample 

2.1. The sample 

We started from the 3D sample of P08. Galaxies were selected 
using J-band absolute magnitudes as a proxy for stellar mass, 
such that Mj < -20.3, which roughly corresp onds to M s te u ar > 
1.5 x 10 10 M o , using a "diet" Salpeter IMF and iBejletal] {2003) 
simplified recipes for deriving stellar mass from J-band lumi- 
nosity. Additional practical constraints were imposed on their 
rest-frame [Oil] equivalent width (i.e., EWo > 15A) and red- 
shift (i.e., z between 0.4 and 0.75), by observing with the multi- 
integral-field unit spectrograph FLAMES/GIRAFFE at the VET. 
This provided a sample of 63 galaxies that represents the J-band 
lumin osity function of galaxies at these redshifts (Yan g et all 
2008). In the following, we restrict our study of the bTFR to 
galaxies within the CDFS (see Sect. 2.4). This sample of 35 
galaxies is still rep resentative of the luminosity function at z~0.6 
(lYang et al.ll2008l) . 

NIR photometry was derived from public images of the 
respective fields where galaxies were selected (i.e., CDFS, 
CFRS, and HDFS). Stellar masses M ste u ar were est i mated from 
M ste [iar I Lk ratios using the method of iBell et all (120031) . as- 
suming a "diet" Salpeter IMF. We estimated in P08 that this 
method, when applied to z~0.6 galaxies, provides us with esti- 
mates of the stellar mass with an associated random uncertainty 
of 0.3 dex, and a systematic uncertainty of -0.2 dex. We note 
that this comparison takes into acco unt the influence o f possible 
secondary bursts of star formation dBorch et al.ll2006l) . In addi- 
tion, the influence of TP-AGB stars on the derivation of stellar 
masses could overestimate the stellar mas s by up to ~0.14 dex 
dMaraston et alj|2006t iPozzetti et alj|2007l) i n a s ystematic way, 
but con stant with redshif t. In lHammer et al.l d2009bl) . it is shown 
that the lBell et al.l d2003) method, when applied to z~0.6 galax- 
ies, provides us with an upper limit to their stellar mass rela- 
tive to the measurements determined by other stellar-population- 
synthesis models and IMF combinations, which is reflected by 
all systematic effects tending to lower these estimates. Using 
these stellar masses to study the TFR is therefore a quite con- 
servative choice, leading to evolutionary trends measured being 
perhaps lower than in reality. We refer the reader to P08 and its 
Appen dix for a clearer descr iption of our adopted method, as 
well as lHammer et al.l (l2009bl) . 



We note that the IMF use d to derive stella r masses (i.e., the 
diet Salpeter IMF following IBell et al.l 12003b differs from the 
IMF used in the next section to derive s tar formation r a tes (i.e ., 
the regular Salpeter IMF). However, in IHammer et al.l (120 09b). 
we argue that stellar masses derived using IBell et al l (120031) 
recipes and a diet Salpeter are roughly equivalent (in a one-to- 
one correlation and not onl y in statistical sense) to ste llar masses 
derived using the models of lBruzual & Charloti (12003) and a reg- 
ular Salpeter IMF. Therefore, in the following, we make no at- 
tempt to correct for this difference in IMF, since the current gen- 
eration of stellar population synthesis models is simply not ac- 
curate enough for difference caused by these two IMFs to be 
significant. 

Compared to P08, we modified the classification of two 
galaxies in the sample. J0332 45. 11-274724.0 was studied in 
detail by IHammer et ail (120091) . who found in HST/ACS im- 
ages from the UDF a dust-enshrouded disk that had not previ- 
ously been seen in GOODS images because of their shallower 
depth (see a comparison between GOODS and UDF images in 
their Fig. 1). Accounting for this disk, the dynamical axis was 
found to be strongly misaligned with the morphological axis. 
IHammer et al.l d2009t) argued that the morpho-kinematic proper- 
ties of J033245. 11-274724. can only be repr oduced by a ma- 
jor merger event. Similarly, Puec h et all (120091) carried out a de- 
tailed study of J033241. 88-274853.9, whose morpho-kinematic 
properties also appeared to be far more closely reproduced 
by a major merger model than that of a simple rotating disk. 
Therefore, we shifted the classification of both objects from RD 
to CK. 



2.2. Star formation rates 

We estimated the total SFR, of of each galaxy, following 
iPuech et alj (l2007bl) . Briefly, SFR to, was taken to be the sum of 
SFRf/y, derived from the 2800A luminosity L2800, and SFR/«, 
which was derived from Spit zer/MIPS photometry at 24yura us- 
ing the ICharv & Elbaz (2001) calibration between rest-frame 
15/zm flux and total IR luminosity Ljr. To co nvert both Ljr an d 
^2800 into SFRs, we used the calibrations of iKennicuttl (1998), 
which rely on a Salpeter IMF. Uncertainties were estimated by 
propagating the flux uncertainty measurement using the SFR cal- 
ibrations. About half of the galaxies in the sample were not de- 
tected by MIPS. In this case, we derived an upper limit to SFR/j ; 
as a function of redshift using Fig. 9 of Le Floc 'h et a El (120051) . 
We assumed that SFR, , is the mean of this limit and SFRj/v/. 
The corresponding uncertainty in this case is half the difference 
between SFRj/y and the limit to SFR/r, which infers a mean rel- 
ative uncertainty of 50% for these objects (see Fig. Q] and Table 

2.3. Gas mass estimates 

To estimate the mass of gas within each galaxy, we used 
the SK relation between SFR and gas densities of IKennicuttl 
(1 19891). which is giv en by S SFS = 2.5 x lO" 4 ^ M Q /yr/kpc 2 . 
iBouche et al.l d2007l) derived a new calibration of the SK law, 
based on distant sources (z~2-3). They found that the SK law 
holds up to z~2.5, but has a quite different power-law index 
(i.e., 1.7 instead of 1.4). For the CDFS sample, we find that 
2g a s ranges from 0.004 to 0.46 M Q /yr/kpc 2 . In this regime, the 
two SK laws are relatively similar, the discrepancy between the 
two relations being caused mostly by the higher star-formation 
densities found in z~2 starbursts (see Fig. 3 of Bouche et al. 
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Fig. 1. SFR as a function of stellar mass in a subsample of 
35 CDFS galaxies. RDs are shown as blue dots, PRs as green 
squares, and CKs as red triangles. The open blue circle corre- 
sponds to the RD+ galaxy for which the velocity measurement 
is more uncertain (see P08). Uncertainties in stellar mass are dis- 
cussed in Sect. 2.1, and are not shown for reasons of clarity. 

12001) . This range corresponds to the star-formati on densities of 
local (U)LIRGs, starburst s, and star- f ormin g disks dBouche et al.l 
120071) . Finally, at z=0.24. lLah et al.l(l2007l) performed direct HI 
measurements of distant galaxies, and found that there is no evo- 
lution in the relation between HI mass and SFR. Therefore, we 
conclude that it is more appropriate to rely on the local SK law 
rather than the one at high-z to derive gas masses in a sample of 
z~0.6 galaxies. 

To do this, we first derived the SFR density "Lsfr = 
S FR tota ilnR 2 gav which was then converted into a gas density 
Y-gas using the SK law. Gas masses were derived as M gas = 
I.g as x nR 2 gas . The same radius R gas , which is defined in the next 
section, is assum ed when normalizing both densities, following 
iKennicuttl (Il989l) . Uncertainties in M gas were derived using stan- 
dard methods of error propagation (see TableQ~]and Fig. It is 
noteworthy that gas masses derived this way are independent of 
the IMF, since the SK law and SFRs are derived using the same 
Salpeter IMjQ. Gas ma sses and gas fractions in z~0.6 galaxies 
are discussed further by Ham mer et al.l(l2009bl) . Of particular in- 
terest here, is that the gas masses derived by inverting the SK law 
lead to a median gas fraction of 31% in z~0.6 galaxies. We note 
that gas fractions, in contrast to gas masses, are not IMF inde- 
pendent. Strikingly, the same value is found using a completely 
different data set and methodology (i.e., based on the evolution 
of the mass-met allicity relation derived from FORS2 long-slit 
spectroscopy, see Rodrigues et al. 2008). It is therefore very un- 
likely that the gas fractions derived at z~0.6 could be affected 

1 The IMF affects the SFR calibration in such a way that both the 
SFR and the numerical factors in the SK relation change and cancel 
eachother. It is possible to derive a direct relati on between UV a nd IR 
luminosities and gas densities, in a similar way to lErb et aT] ||2006) using 
the Ha luminosity, therefore eliminating the SFR density, which is the 
IMF-dependant quantity. 



Fig. 2. Gas mass as a function of stellar mass in a subsample of 
35 CDFS galaxies. RDs are shown as blue dots, PRs as green 
squares, and CKs as red triangles. The open blue circle corre- 
sponds to the RD+ galaxy for which the velocity measurement 
is more uncertain (see P08). Uncertainties in stellar mass are dis- 
cussed in Sect. 2.1, and are not shown for ease of clarity. The 
blue line represents equality between the stellar and gas masses, 
while the blue dashed line re presents the local re lation between 
these two quantities given by Schiminovich (2008). The median 
gas fraction is found to be 31%, or M gas /M sle n ar =0A5, which is 
represented as a thin black line. 



by any systematic effect, which makes us quite confident in our 
estimates of gas masses. 

2.4. [Oil] gas radius estimates 

To convert from density to mass, one needs to define and mea- 
sure a gas radius R gas . To do this, we first constructed rest-frame 
UV images by summing observed B and V bands. Given that 
not all galaxies in the sample have homogeneous imaging (see 
iNeichel et alj2008l) . we decided to restrict the study of the bTFR 
to galaxies with the highest quality images. Therefore, we re- 
stricted the study of the bTFR to the subsample of 35 galaxies 
lying in the CDFS with HST/ACS images, which is still repre - 
sentative of the luminosity function at z~0.6 dYang et al.f 2008). 
This allowed us to limit uncertainties in the derivation of gaseous 
radii and therefore, gas masses. 

From rest-frame UV images, we derive d for each galaxy 
an ax is ratio b/a and a PA using SExtractor dBertin & Arnoutsl 
1996). For each galaxy, these two parameters were used to gen- 
erate a set of images with flat ellipses of increasing radii using 
the IDL procedure DIST_ELLIPSE. To account for the relatively 
coarse spatial sampling of the GIRAFFE IFU (0.52 arcsec/pix, 
which roughly corresponds to 3.5 kpc at z~0.6), these high- 
resolution ellipses were then rebinned to the GIRAFFE pixel 
scale after determining as accurately as possible the position 
of the GIRAFFE IFU on the UV image (see, e.g.. iPuech et al.l 
l2007bl) . These simulated [Oil] GIRAFFE images were renor- 
malized in terms of flux using the observed IFU-integrated [Oil] 
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value, and pixels with a resulting flux lower than the mini- 
mal [Oil] flux detected within the GIRAFFE IFU for a given 
galaxy were disregarded, to account for GIRAFFE IFU pixels 
being selected in terms of signal-to-noise ratio (i.e., S/N>3, see 
Flor es et al.ll2006l) . By identifying the simulated [Oil] map that 
most closely reproduces the observed one, it is in principle pos- 
sible to retrieve the underlying [Oil] total-light radius. However, 
the exact position of the IFU, as well as the precise value of the 
seeing during the observa tions are subject to uncertainties (e.g., 
iPuech et alj|2007bl 120081) that prevent us fr om being able to de- 
termine the exact position of the IFU grid. The relatively coarse 
spatial resolution of the GIRAFFE IFU also makes it relatively 
difficult to retrieve precisely the underlying [Oil] flux distribu- 
tion. To mitigate the uncertainty associated with these effects, 
we did not use a classical chi-square minimization between the 
simulated and observed [Oil] maps, but used a spatial cross- 
correlation to determine, for each galaxy, the subset of simulated 
[Oil] maps that maximizes the number of pixels illuminated in 
both the observed and simulated [Oil] maps. This subset of sim- 
ulated maps corresponds to a range of ellipse radii that are de- 
convolved from the GIRAFFE coarse spatial sampling. We as- 
sumed that the mean value of this range is a measure of the total 
[Oil] radius R[pir\> whose error-bars are equal to half this range. 
Finally, R[on] were quadratically deconvolved from a 0.8 arcsec 
mean seeing disk. The maximal value of the ratio of the error- 
bars to the corresponding radii is found to be 51%, while their 
distribution shows an average of 17% and a l-cr spread of 12%. 
The derived radii and errorbars in kpc are shown in Fig. |4] and 
listed in Table Q] 

3. Evidence that major mergers are responsible for 
the large scatter in the distant TFR 

3.1. The z~0.6 stellar-mass TFR 

The smTFR in the IMAGES sample was initially derived in P08 
(see their Appendix). We show in Fig. [3] (see left panel) a re- 
vised version of this relation, in which the classification of two 
objects has changed (see Sect. 2). We also show here all galaxies, 
including non-relaxed ones. The dispersion in this relation was 
shown to be caused entirely by non-relaxed systems by P08 (see 
next section). Maintaining the slope at its local value, we found a 
shift in the smTFR zeropoint of 0.34^ ^ dex in stellar mass, or, 
equivalently, of ~0. 12 dex in velocity, between z ~ 0.6 and z = 0, 
consistent with the shift derived by P08. The errorbars account 
for possible systematic effects, which were quantified in detail 
by P08. Most of them are associated with the estimation of stel- 
lar mass. As we noted in Sect. 2.1, all these effects would tend 
to decrease the estimates of the stellar mass i n z~0.6 galaxies , 
since we "maximized" them by adopting the iBell et all d2003l) 
method. The only systematic effects that could reduce the shift 
in zeropoint are relatively insignificant (i.e., -0.06 dex), and are 
associated with the derivation of rotation velocities (see P08). In 
P08, we found that the shift of the smTFR with redshift cannot 
be interpreted as a pure evolution along the velocity axis, and 
that only half, at most, could be accounted for by a velocity shift 
caused by gas accretion within the optical radius. We then con- 
cluded that most of the shift in the smTFR reflects an increase in 
ste llar mass in rotating disks over the past 6 Gyr. 

iKoen & Lombard! d2009l) performed a permutation test to 
demonstrate that the distant K-band TFR found by P08 was off- 
set significantly from the local relation. We repeated their anal- 
ysis for the smTFR and found a probability <k1% that the two 
relations have identical slope and zeropoint. To test the reliability 



of this result in terms of systematic uncertainties, we shifted the 
distant smTFR in steps of 0.01 dex in stellar mass toward the lo- 
cal relation and repeated the permutation test. We found that for 
this probability to exceed 10% we would need to bias the stellar 
mass in the distant sample by +0.28 dex, which exceeds system- 
atic uncertainties (see above). This means that, within the sys- 
tematic effects on the zeropoint shift identified by P08, the two 
relations differ significantly, which secures the results of P08. 

3.2. Origin of the scatter in the z~0.6 smTFR 

In iFlores et"aTj d2006l) . we showed that the dispersion in the 
z~0.6 smTFR was caused by non-relaxed galaxies. This result 
was confirmed by P08, who showed that the dispersion in the 
distant smTFR, once restricted to RDs (i.e., 0.12 dex), was sim- 
ilar to that in the local relation, with 0.15 dex. The galaxies con- 
tributing to the scatter in the distant s mTFR exhibit stro ng per- 
turbations in their kinematic maps (e.g jYang et al.l2008l) . which 
led us to claim that most of them were probably associated with 
major mergers (see P08). This interpretation is also supported 
by the large scatter seen in the specific angular momentum ver- 
sus rotation velocity plane , which cou l d be produced by these 
events, as demonstrated by P uech et alj d2007l) . 

Interestingly, IWeiner et alj d200 6) defined a new kinematic 

tracer S = Jo.5 x V^ ot + & 1 that combines rotation velocity and 
velocity dispersion. They found that velocity dispersion is an 
important component of the dynamical support in intermediate- 
redshift galaxies, especially for morphologically compact or dis- 
turbed systems. T heir study confirmed the earlier finding of 
IPuech et ail ((2006) that the majority of compact galaxies have 
motions in which dispersion plays an important role (see also 
IPuech et alj2007l) . lKassin et al.ld2007l) later showed that using S 
in the smTFR leads to a signifi cant reduction in its scatter. In ad- 
dition, Covington et al. (2009) showed, using hydro-dynamical 
simulations of major mergers, that this reduction in scatter is 
consistent with a (partial) transfer of energy between the kinetic 
energy associated with bulk motions to that associated with ran- 
dom motions. This transfer is driven by s hocks and collisions 
gener ated during the merger events (e.g., IMontero-Ibero et ail 
l2006h . 

That the "merger-hypothesis" is the underlying cause of the 
large scatter measured in the distant TFR can be tested by search- 
ing for specific objects and verifying whether this transfer of 
energy between bulk and ra ndom motion s drive n by interaction- 
driven shocks really occurs . IPeirani et aT] (l2009) studied one CK 
galaxy of the sample (J033239.72-275 154.7) in detail and they 
demonstrated that both its morphology and kinematics were re- 
produced well by a 1:3 major merger where the companion 
galaxy had been in a retrograde parabolic orbit and there had 
been an inclination angle of 15 degrees between the two or- 
bital planes. This object is a perfect test-bed for the present dis- 
cussion, since it is measured to have log V// fl ,= 1.48 km/s and 
M stellar — 10.31 M , i.e., to be among the CK galaxies in the dis- 
tant smTFR that ar e the most deficient in rotation velocity (see 
left panel of Fig. [3}. Peir ani et al.l |2009) showed that, because of 
the retrograde nature of interaction, the host galaxy can lose an- 
gular momentum to its companion, resulting in the deceleration 
of the main progenitor. This explains why the resulting merger 
remnant appears to have systematically smaller velocities in the 
smTFR than other objects. Shocks in the gaseous phase gener- 
ated during the interaction produce an increase in the velocity 
dispersion, similar to that suggested by Coving ton et alj (2009) 
by their merger simulations (see their Fig. 1). If the scatter in the 



Table 1. Principle properties of the sample of 3 5 galaxies used in this study, in order of increasing RA (see Sect. 2.5). From left to right: IAU ID, redshift (from Ravikum ar et al.l 
|2007[) . dynamical class (from [Yang et al.ll2.Q08l) . rotation velocity and associated uncertainty (from P08), stellar mass log(M sle u ar /M e ) (from P08), SFRyy (in M G /yr), SFR /fl 
(in M /yr), total SFR and associated uncertainty ASFR (in M /yr), rest-frame UV half-light radius Ryy and associated uncertainty ARf/y (in kpc), [Oil] total radius R[o//j and 
associated uncertainty AR[o//j (in kpc), and mass of gas log (M gas /M Q ) with associated uncertainty A log (M ste n ar /M Q ). 
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Fig. 3. left: Evolution of the stellar-mass TFR in the sample of 64 galaxies of P08 (RDs are shown as blue dots, PRs as green 
squares, and CKs as red triangles). The open blue circle correspond to the RD+ galaxies for which the velocity measurement is 
more uncertain, see P08). The two open red triangles correspond to the two galaxies that were changed from RD to CK compared 
to P08 (see Sect. 2.1). The black line is the local smTFR, while the blue dash-line represent a linear fit to the z ~ 0.6 smTFR. Left: 
Same relation but using S . The black and blue-dashed lines are the same as in the left panel. 



smTFR is really caused by a transfer of energy between bulk and 
random motions driven by merger events, then this object should 
fall back onto the TFR, once its velocity dispersion is taken into 
account. 

To test this, we derived S in the IMAGES sample u sing 
Vrot - V 7 fiat and <x = ov,^, as defined by |Puech~et al. (2007). We 
show the resulting 5 -TFR in Fig. [3](see right panel). Using S re- 
sults in the shifting of all RDs back onto the local smTFR, so we 
propose to use the local slope in all future investigations. We also 
note that several CK galaxies lie at relatively high S compared to 
the bulk of the sample. These galaxies were found to have large 
rotation velocity in the conventional smT FR. However , these 
galaxies do not s eem to be present in th e [Kassin et"ak I (12007b 
smTFR, although Covington et al] d2009t) managed to produce 
some of these galaxies in their merger simulations (see their 
Fig. 3). These differences seem to be related to different selec- 
tion criteria involving inclination. We did not apply any selec- 
tion of this type (see P08), while Kassin et al. (2007) removed 
galaxies with low inclinations (i.e., inc < 30deg) because of 
the larger associated uncertainty, and those with high inclina- 
tions (i.e., inc ^ 70de g), because of the effect of dust (see also 
ICovington etai]|2009l) . 

The total scatter in the resulting S -TFR is significantly lower 
than in the smTFR, since it decreases from 0.63 to 0.34 dex. This 
residual scatter is still a factor of two larger than the residual dis- 
persion in the local relation, which is cr rcs =0.15 dex (in mass). 
Looking at the dispersion in more detail, the scatter remains 
roughly constant in RDs and PRs (from 0.12 to 0.14 dex and 
0.32 to 0.29 dex, respectively), but strongly decreases from 0.83 
to 0.41 dex among CKs. That kinematically disturbed galaxies 
show a smaller but still large scatter compared to the local re- 
lation suggests that part of the remaining discrepancy is related 
to the greater uncertainty in rotation velocity of CKs (see P08). 
The median uncertainty in the rotation velocity of the sample is 



found to be 0.12 dex, which translates into ~0.3 dex in stellar 
mass, i.e., the same order of magnitude as the residual scatter in 
the distant S -TFR. 

As one can see in the right panel of Fig. [3] accounting for 
velocity dispersion brings the merger-test J033239. 72-275 154.7 
back onto the S -TFR. This is the first direct o bservational ev- 
idence of the transfer of energy described by ICovington et alj 
(2009). As shown above, this process can account for most of the 
scatter in the smTFR. Given the still relatively large uncertain- 
ties, we cannot exclude other process(es) playing a role in pro- 
ducing the scatter in the distant relation, but these results clearly 
indicate that this transfer of energy is the main underlying cause, 
and that mergers are the main cause o f the large sc a tter se en in 
the z~0.6 TFR, as initially argu ed by iFlores et"aT] 12006) (see 
also iKannappan & Bartonll2004l) . 

4. Gas content of intermediate-mass galaxies 

4. 1 . Comparison between gas and stellar extents 

We show the total-light [Oil] radius R[oii] as a function of 1.9 
times the UV half-light radius Ruv in Fig. [4] For a thin expo- 
nential disk, 1.9 time s the half-light radius is equal to the optical 
radius of the galaxy dPersic & Salucci|[l991l) . which provides us 
with a useful comparison poinfl Galaxies that have a R[oii] sig- 
nificantly lower than 1.9xRuv are galaxies for which the UV 
light extents farther than the IFU FoV (3x2 arcsec 2 ). Two ex- 
amples of these galaxies can be seen on the right side of Fig. [4] 
For these galaxies, it is clear that assuming R[oii] to be a mea- 
sure of the gaseous extent would lead to underestimating the gas 
mass. Therefore, we defined R gas = MAX(R[on], 1-9 x Ruv) to 

2 For a thin exponential disk, this radius equals 3.2 times the disk 
scale-length, which defines the total stellar-li ght radius, and is statisti - 
cally equivalent to the isophotal radius R 2 5, see lPersic & Saluccli 19911) . 
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Fig. 4. Comparison between the [Oil] and UV radii for the CDFS 
sub-sample of galaxies. The two insets on the left (30-arcsec 
wide) show MIPS imaging at 24 fim of J033228.48-274826.6 
(upper one) and J033244.20-274733.5, while the two insets 
on the right show z-band HST/ACS images of J033230.78- 
275455.0 and J033226.23-274222.8 superimposed with the 
GIRAFFE IFU grid. RDs are shown as blue dots, PRs as green 
squares, and CKs as red triangles. The open blue circle corre- 
sponds to the RD+ galaxy for which the velocity measurement 
is more uncertain (see P08). 



be a conservative measure of the gaseous extent. We note that 
for galaxies with non-relaxed kinematics, this might still under- 
estimate the gas radius. The uncertainties correspond to the re- 
spective uncertainties in R[oii] or 1.9xRuv depending on which 
radius was the largest for a given galaxy. 

Most RDs (see blue dots) fall relatively close to the equality 
line, which instills confidence in the estimate of R gas in dynami- 
cally relaxed systems. To first order, one can assume that, at least 
on large spatial scales and for relaxed systems, the spatial distri- 
bution of UV light emitted by young OB stars is correlated with 
the emission from the [Oil] ionized gas detected by GIRAFFE. 
However, in more disturbed galaxies, there is a clear trend in 
which the gas is more extended than the UV light. We find a 
median [OII]-to-UV radius ratio of 1.3+0.16 (1-cr bootstrapped 
uncertainty), which is consistent with the mean r atio of ~ 1 .2 of 
emiss ion line to B-band scale-lengths found by Ba mford et alj 
(2007) in a sample of field galaxies at similar redshifts. Using a 
Student t-test, the probability that this median ratio is equal to 
one (i.e., that the total [Oil] and UV radii are statistically equiv- 
alent) is found to be 6%. This trend is even more pronounced 
for more compact systems: galaxies with a UV radius lower or 
equal to the median in the sample have a median ratio of ~1.7, 
while those having a UV radius larger than the median value, 
have a ratio -0.9. We also examined the GOODS/MIPS (DR3) 
images of the four galaxies with a UV radius smaller that one 
GIRAFFE pixel. All were detected, although only three at a sig- 
nificantly high level for their flux to be measurable. Interestingly, 
two of them were resolved by MIPS (see the two insets on the 
left side of Fig. [U, which confirms that the gas content can be 
significantly more extended than the UV stellar light. 



4.2. What causes the large gaseous extent in z~0.6 
galaxies? 

It is interesting to investigate what causes this discrepancy be- 
tween the extents of the ionized gas and the UV ionizing stellar 
light from OB stars. Most of the ionizin g stellar light is as sumed 
to be emitted by stars earlier than B2 ( Strombergl l 1 9391) . so to 
first order, one would expect conversely that the ionized gas is 
confined to within the UV radius. What mechanisms can explain 
our actual findings? 

Supernovae (SN) and AGN feedback are two mechanisms 
that are known to be able to drive outflows in galaxies. If these 
mechanisms were responsible for ejecting ionized gas to the UV 
radius, we should detect offsets between absorption and emis- 
sion lines in a large fraction of galaxies in the sample. For 
20 galaxies in the GIRAFFE s ample, we were able to retrieve 
FORS2 integrated spectra (see iRodrigues et al.l 120081) that al- 
lowed us to measure absorption lines. Among these galaxies, 
we found systematic shifts in only three or possibly four of 
them(ata~100km/slevelinJ033224.60-274428.1,J033225.26- 
274524.0, J033214.97-275005.5, and possibly in J033210.76- 
274234.6). Hence, we conclude that outflows, regardless of their 
powering mechanism, cannot explain why gas extends farther 
out than the UV stellar light. 

If ionized gas is not pushed out of the UV disk by outflows, 
then the gas must be ionized by a mechanism other than the ra- 
diation of OB star s. Interestingly, for one object in the sample, 
iPuech etail J2009) observed a region that is completely devoid 
of stars, where ionized gas was detected. They concluded that 
this gas was ionized by shocks induced by a major merger. These 
interaction-driven shocks c ould be an important m ean of ioniz- 
ing gas in z~0.6 galaxies. iHammer et alj d2009bl) showed that 
the morpho-kinematics of most galaxies in the sample can be re- 
produced by the simulations of major mergers. We note that gas 
in all galaxies in the sample have lower V/cr ratios than their lo- 
cal cou nterparts, which m ight be a dynamical signature of these 
shocks dPuech et al.l 120071) . In partic ular, most PRs show off- 
center velocity dispersion peaks (see Yan g et al]|2008l) . which 
strengthens this interpretation. 



5. The baryonic Tully-Fisher relation at z~0.6 

5.1. The local relation 

We used the local sample gat hered by IMcGaughl d2005l) to de- 
rive our local reference bTFR. iMcGau gh (2005) u 



(amongst others) the iBell et al" 



McGaugh (2005) used V flat and 



(2003) method to estimate stel 
lar mass, so their sample is particularly well suited to compare 
with our high-z bTFR. However, they used the B-band lumi- 
nosity to derive stellar mass, which is a less accurate tracer 
than K-band luminosities. Hence, we rederive d stellar ma s s from 
K-band luminosities by cross-correlating the IMcGaughl {2005) 
sample with the 2MASS database. Stellar masses were derived 
following the method outlined in P08 for their local sample, i.e., 
by accounting for k-corrections and extinction. We found that 
log (M baryonic /M Q ) = 2.10 + 0.42 + (3.74 + 0.20) x log (V flal ) has 
a residual scatter of cr= 0.25 dex. In the f ollowing, we checked 
that by instead using the lMcGaughl (f2005) stellar masses derived 
from the B-band luminosity (see his Table 2 with P-l) would 
not c hange significan tly our results. 

IStark et all d2009h established a new calibration of the lo- 
cal bTFR using gas-dominated galaxies, which has the advan- 
tage of being l ess sensitive to a given estimator of the stellar 
mass. However. lStark et ail ([2009) averaged the bTFR produced 
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by various stellar mass estimators. Therefore, we chose to keep 
our own derived relation, which has the advantage of resulting 
from the same methodology used to derive the smTFR, both at 
low and high redshift, as well as the bTFR at high redshift, which 
minimizes possible systematic effects. We nevertheless note that 
our slope and zeropoint are consistent with their calibration. In 
contrast to the smTFR, whose calibration depends on the galaxy 
populations analyzed, the bT FR is found to be independent of 
galaxy type (IStark et al]|2009l) . We therefore do not expect any 
bias to be caused by the choice of a particular local sample and 
its level of representation in the local population. 



5.2. The baryonic TFR at z~0.6 

We show the first bTFR obtained at high redshift so far, in Fig. 
[5] Overall, the bTFR has the same structure as the smTFR, with 
non-relaxed galaxies associated with the greatest amount of dis- 
persion. Holding the slope constant, we found a distant zeropoint 
of 2.16 + 0.07, i.e., a shift of +0.06 dex in baryonic mass be- 
tween the z=0 bTFR and z~0.6 RDs. This shift translates into 
a -0.02 dex shift in rotation velocity. The scatter in the distant 
bTFR for RDs is found to be roughly similar to that in the local 
bTFR, with 0.14 dex, which suggests no evolution in slope. In 
comparison with the distant smTFR, the scatter in RDs is found 
to be roughly equivalent, with 0.14 dex in the bTFR instead of 
0.12 dex in the smTFR (i.e., bootstrapping reveals that the l-cr 
spread of the scatter in RDs in the smTFR is ~0.03 dex). In con- 
trast, the total scatter is found to be significantly larger, with 0.83 
dex in the bTFR instead of 0.63 in the smTFR. We used Monte- 
Carlo simulations to investigate whether the larger uncertainty 
associated with the gas mass, compared to the stellar mass, is 
responsible for the increase in scatter in the bTFR. To do this, 
we simulated 1000 smTFRs by moving the galaxies within the 
errorbars in their baryonic mass. We then refitted the resulting 
smTFRs, and found that the total scatter (as well as the scatter 
for RDs) does not increase significantly. Therefore, it appears 
that the bTFR has a significantly larger scatter than the smTFR, 
at least for the total galaxy population. 

In assessing whether there is an evolution in the intercept of 
the bTFR, it is important to account for possible systematic ef- 
fects, which are usually more significant than random effects in 
the TFR (see P08 and Sect. 3). In the following, we express all 
systematic uncertainties in terms of their influence on the evo- 
lution of the zeropoint of the bTFR between z~0.6 and z=0. 
Systematic effects on the velocity estimates in the distant sam- 
ple were found to be +0.02 dex, which corresponds to ±0.08 
dex once converted into Mi, myo „i C . Regarding the stellar mass, 
we identified in P08 two possible systematic effects, which are 
an evolution in the IMF (+0.05 dex), and an evolution in the K- 
band mass-to-light ratio (+0.1 dex) with redshift (see Appendix 
of P08). Finally, concerning gas masses, the only possible sys- 
tematic effect could be an evolution of the SK law with redshift, 
since we argued in Sect. 2.3 that this quantity is derived inde- 
pendently of the IMF. As stated in Sect. 2.3, there is no evidence 
of an evolution in the SK law between z~0.6 and z=0, based on 
a comparison between their star-formation densities and those 
of local starbursts. Therefore, the main systematic effect that 
could affect the baryonic mass would be that related to stellar 
mass, which translates into a possible impact on the evolution of 
the bTFR zeropoint shift between z~0.6 and z=0 by +0.15 dex 
(0.05+0.1, see above). In conclusion, we find a bTFR zeropoint 
shift between z~0.6 and z=0 of -0.05+0.08 (random, see above) 
-0 08 ( s y stema ti c )- 
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Fig. 5. Evolution of the baryonic TFR in the C DFS subsampl e 
(see text). The black line is the local bTFR of lMcGaughl d2005l) . 
while the blue-dashed line is the fit to the distant relation. All 
panels have similar limits to ease the comparison between the 
different relations 



Using the permutation test of iKoen & Lombard! d2009l) . we 
find a probability of ~62% that the local and distant relations 
have the same slope and intercept. Given that there are only 
six RDs in the distant bTFR, we also tested whether the local 
and distant intercepts are different, assuming that the slope does 
not ev olve. To do this, we used a Welch t-test (IKoen & Lombard! 
2009), and found that the probability that they are equal is 91%. 
We also found that this probability drops below 10% when we 
systematically increase the local bTFR zeropoint by at least 0.77 
dex: this means that the local and distant bTFR zeropoints are not 
statistically different within systematic uncertainties (see above), 
assuming no evolution in slope. As stated above, the residual 
scatters of the local and distant relation are similar, which sug- 
gests that assuming there is no evolution in slope is a reasonable 
hypothesis. Therefore, we conclude that we do not detect any 
significant shift in the bTFR between z~0.6 and z=0. 

6. Discussion & conclusions 

There is now little doubt that the large scatter found in the z~0.6 
TFR is cause d by non-relaxed systems associated with ma- 
jor m e rgers dFlores et alj 20061 : iKassin et al.1 120071: iPuech et al" 
l2Q08h ICovington etaiT T2009: see also iKannappan & Barton! 
|2004|) . Mergers provide a natural and coherent frame for inter- 
preting the morpho-ki nematics properties of z~0.6 galaxies (see 
iHammer et alj|2009bl) . In this paper, we have reported that ion- 
ized gas extends farther out than the ionizing light from OB stars: 
we indeed found a probability of 6% that the median ratio of the 
gas to stellar radius is equal to one. This can be seen as a nat- 
ural consequence of shocks produced during these events. It is 
even difficult to figure out how gas lying significantly farther out 
than any radiation source could be ionized by any other process, 
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given that a very negligible fraction of objects in the sample are 
found to contain evidence of gas outflows, which might spread 
ionized gas out to the UV radius. 

Returning to the TFR, we confirm that, once restricted to 
rotating disks, the z~0.6 smTFR appears to be shifted toward 
fainter masses by 0.34 dex. This shift is found to be very sig- 
nificant within the random and systematic uncertainties. In con- 
trast, by including the gas fraction of the baryonic mass in the 
TFR, we found that, within the uncertainties, the z~ 0.6 relation 
is con sistent with the local relation derived for the McGaugh 
(2005) sample. This implies that star formation in rotating disks 
is mostly fed by gas that is already gravitationally bound to 
galaxies, otherwise we would have detected significant evolu- 
tion in the bTFR. This does not mean that there is no external 
gas accretion: given the still relatively large associated uncer- 
tainties, there could be room for external gas accretion, at a level 
of up to roughly one third of the local baryonic mass. However, 
there is presently no nee d for this ext ernal gas accretion. 

At higher redshifts, ICresci et"aT | (120091) found an evolution 
of 0.41 dex in the zeropoint of the smTFR between z~2.2 and 
z=0, i.e., of similar amplitude (within uncertainties) to the evo- 
lution found in this paper between z~0.6 and z=0. This would 
imply that the smTFR does not evolve significantly between 
z~2.2 and z~0.6, while it evolves significantly between z~0.6 
and z=0. To try to unde rstand this surprisi ng result, we refit- 
ted the smTFR found by ICresci et~aT] (120091) . but assumed that 
the local slope of reference is that used in the present study. 
The choice of the slope in the local relation is crucial for de- 
riving the evolution in zeropoint as a function of redshift (see 
Sect. 5.1 of P08 fo r details). While the distant sample used by 
ICresci et~ai] (|2009) seems to be representative of z~2 galaxies, 
at least for galaxies with stella r masses higher than ~2xlO lo M 
dForster Schreiber et aljl2009l). they used the local relation de- 
rived by iBell & de Jongl d200ll) . which relies on the sample of 
IVerheiienl d2001l) that was shown to be biased toward an excess 
of low-mass, gas-rich galaxies by Ha mmer et al.l d2007l) . In P08, 
we derived a lo cal smTFR using a re presentative subsample of 
the SDSS from Pizagno et al] (120071) . In this sample, the slope 
is found to be smaller, i.e. , 2.8 (see Appendix of P08) instead 
of the val ue of 4.5 used by Cresc i~et~aT] d2009l) . By refitting the 
relation of ICresci et alj d2009l) using the slope derived from this 
representative local sample, we find an evolution in zeropoint of 
~0.6 d ex, instead of ~0.4 dex found using the Bell & de Jong ] 
2001) slope. We conclude that the shift found by ICresci et al.1 



2009) probably underestimates the evolution in zeropoint be- 
tween z~2. 2 and z=0. 

Finally. ICresci et al.l (120091) interpreted the shift of the z~2.2 
smTFR as the result of gas accretion onto the forming disks in 
filaments and cooling flows, as suggested by theoretical and nu- 
merical models. This gas accretion process could also play a role 
in feeding the outer regions of z~0.6 galaxies with fresh gas, to 
within a limit of 30% of the local baryonic mass, as discussed 
above. However, it would remain unclear what the ionization 
source of this gas is (see Sect. 4). These models also clearly pre- 
dict that these cold flows are strongly suppressed in the massive, 
z~0.6 haloes th at the galaxies studied in this paper inhabit (e.g ., 
iDekel et al.l2009t) . In a companion paper (iHammer et al. 2009b), 
we explore another possibility that a large fraction of local spi- 
rals could have rebuilt their disk following a major merger at 
z<l. This scenari o, dubbed as "spiral r ebuilding disk scenario", 
was proposed by Ham mer et al.l d2005l) . owing to the remark- 
able coincidence of the evolution of the merger rate, morphol- 
ogy, and fraction of actively star-forming galaxies. In this sce- 
nario, major mergers expel gas in, e.g., tidal tails, which is later 



re-accreted to rebuild a new disk ((Barnes 2Q02 iRobertson et al.l 
120061 iHopkins et al.l 12009b . It is tempting to associate a sig- 
nificant part of the gas reservoir with this process, which can 
therefore provide us with an evolutionary framework in which 
to intepret the non-evolution of the bTFR with redshift, and the 
finding that ionized gas extends farther out than the stellar UV 
light. A dynamical imprint of this accretion might already have 
been detected in terms of the lowe r V/cr found in distant gaseous 
disks compared to local spirals dPuech et al.ll2007l) . Therefore, 
different physical mechanisms could be driving galaxy evolution 
and the shift in the smTFR, depending on redshift, which would 
complicate the interpretation of the evolution in the smTFR from 
very high redshift to z=0. 
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